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ABSTRACT 

The extremely luminous supernova 2006gy (SN 2006gy) is among the most energetic ever 
observed. The peak brightness was 100 times that of a typical supernova and it spent an un- 
heard of 250 days at magnitude -19 or brighter. Efforts to describe SN 2006gy have pushed 
the boundaries of current supernova theory. In this work we aspire to simultaneously repro- 
duce the photometric and spectroscopic observations of SN 2006gy using a quark nova model. 
This analysis considers the supernova explosion of a massive star followed days later by the 
quark nova detonation of a neutron star. We lay out a detailed model of the interaction be- 
tween the supernova envelope and the quark nova ejecta paying special attention to a mixing 
region which forms at the inner edge of the supernova envelope. This model is then fit to pho- 
tometric and spectroscopic observations of SN 2006gy. This QN model naturally describes 
several features of SN 2006gy including the late stage light curve plateau, the broad Ha line 
and the peculiar blue Ha absorption. We find that a progenitor mass between 20M Q and 4OM 
provides ample energy to power SN 2006gy in the context of a QN. 

Key words: Dense matter - Radiative transfer - Stars: evolution - (Stars:) supernovae: indi- 
vidual: 2006gy 



1 INTRODUCTION 

The supernova (SN) 2006gy discovered by Robert Quimby in Au- 
gust 2006 has challenged our understanding of stellar evolution 
( |Quimby et al^|2007) . SN 2006gy was 100 times more luminous 
than a typical SN and at the time the most energetic ever recorded. 
For almost a year it continued to radiate at a pace in which an ordi- 
nary SN could only sustain for at most a few days. Explaining this 
tremendous energy budget pushes the limits of current SN theory. 

SN 2006gy reveals many singular features in both its light 
curve (left panel of Fig. [TJ and Hydrogen spectrum (right panel 
of Fig. [TJ. The light curve exhibits a luminous peak, broad shape 
and energetic plateau while the Hydrogen spectrum displays a curi- 
ous evolution and a peculiar structure. Models have been proposed 
to describe individual characteristics of this event but all have been 
left wanting. 



1.1 Supernova Explosion Models 

The significance of SN 2006gy was first discussed independently 
by 
by 



Smith et al.|<|2007fr and |Ofek et"aL| l |2007[ >. The model proposed 



Ofek et al. ( |2007 1 involved a type la SN exploding during the 



common envelope phase of a binary system. The collision between 
the SN ejecta (SNE) and dense circum-stellar material (CSM) re- 
leases the amount of energy required to explain the observations of 
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SN 2006gy. However this model demands the CSM to be extraordi- 
narily massive (Ofek et al. 2007 1 suggesting a mass loss rate several 
orders of magnitude greater than expected ( |Yungelson et al.|2008| >. 
In addition the spectrum of SN 2006gy indicates the presence of 
elements not seen in a type la SN ( |Smith et al.|2010| . 

|Smith et al.| ( |2010| > championed a CSM model which consid- 
ered a wind from an exotic luminous blue variable (LBV) to ac- 
count for the massive CSM. Although this model can account for 
the first five months of the observed light curve, the diffusion pro- 
cess requires a rapid drop off in luminosity ( Agnoletto et al.|2009| > 
rather than the observed plateau. Attempts to reconcile the CSM 
model with the late stage plateau invoke the decay of radioactive 
56 Ni and 56 Co to generate the needed luminosity. This necessitates 
10 to 27 times the maximum amount of 56 Ni that can be created 
by a SN l jUmeda &"N omoto 2008). This description for the late 
stage light curve of SN 2006gy has been rebuked by recent near- 
infrared observations which s how that the decline in luminosity is 
inconsistent with 56 Co decay ( |Miller et al.|2010| >. A further chal- 
lenge to the CSM model lies in the fact that a high mass-loss rate is 
deduced using an equation that requires the stellar wind to be con- 
stant in time; however the conclusion made by |Smith et al.|{2010[ > 
is that the stellar wind must vary with time, contradicting their ini- 
tial assumptions (Dwarkadas 2011 1. For this model to achieve the 
amount of radiated energy observed in SN 2006gy, the initial ki- 
netic energy of the SNE must be at least 5 x 10 5 ' erg. This implies 
that SN 2006gy was one of the most energetic SNe and requires a 
massive (M CSM > 20M S ) CSM cloud ( |Smith et al.|2010fr . As noted 
by |Smith et al.|(20T0[ >, the transformation of kinetic energy into ra- 
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diation should translate into substantial narrowing of the Ha line 
as the fast moving ejecta slows, which is contrary to that observed. 
Finally any model that involves the collision of SNE with a dense 
CSM should be a strong emitter of X-ray radiation since the shock 
temperature would be very high (Blinnikov 2008). Observations 
from the Chandra X-ray observatory of SN 2006gy have seen very 
limited X-ray emission ( |Smith et al.p OlOl. The CSM model ex- 
plains the lack of X-rays through self-absorption by the cold, outer 
layers of the CSM ( |Blinnikov|2008l >. 

A pulsational pair-instability (pPISN) model for SN 2006gy 
was considered by Woosley et al. (2007). In this model an unusu- 
ally massive star ( > 100 M Q ) becomes prone to the y = 4/3 insta- 
bility, triggering a SN explosion. The pair-instability process must 
occur twice, leading to a collision between ejected shells which re- 
leases the energy needed for SN 2006gy. In order to reach the peak 
luminosity, it was necessary for Woosley et al. ( 2007 ) to artificially 
increase the kinetic energy of the second ejection. Like the CSM 
model, the pPISN model fails to achieve on several levels. First it 
is unable to properly account for the light curve plateau, as it falls 
off too rapidly ( Woosley et al. 2007 1. Woosley et al. ( 2007 1 explain 
that the extra energy required to fit the plateau must by generated by 
radioactive decay of 56 Co, a conclusion contradicted by late stage 
near-infrared observations ( Miller et al. 2010 1. Secondly, the multi- 
component structure of the Ha line is difficult to reconcile with the 
pPISN model which demands all the Hydrogen to be contained in 
an outer shell ( Woosley et al. 2007 ; Blinnikov 2008 ). Thirdly, a pro- 
genitor star of 1 10 M is required for the pPISN model to achieve 
the power output of SN 2006gy. A star this massive is expected to 
lose its Hydrogen long before it goes SN (Yung elson et al.|2008| >. 
In order for numerical simulations to create a massive star with 
sufficient Hydrogen the merger of several young high mass stars 
was required (Yungelson et al.|20 08l. This type of run-away stel- 
lar merger would be common in the early Universe but atypical 
240 million years ago ( Yungelson et al. 2008 1, at the era of the SN 
2006gy explosion. 

In this paper we present the quark nova as an alternative ex- 
planation for SN 2006gy. Section 2 describes the physics of the 
interaction between the supernova envelope and quark nova ejecta. 
In section 3, the photometric and spectroscopic fits of the quark 
nova model to observations of SN 2006gy are presented. Section 
4 provides a discussion of other observations of SN 2006gy which 
are consistent with the quark nova model. Finally a conclusion is 
presented section 5 and future work discussed in section 6. 



2 QUARK NOVA MODEL 

The quark nova (QN) was proposed as an alternative explanation 
for SN 2006gy ( |Leahy & Ouyedl|2008"l |Ouyed et al!]|2009a) . A 
QN is expected to occur when the core density of a neutron star 
reaches the quark de-confinement density and triggers a violent 
( |Ouyed et al.||2002) conversion to the more stable strange quark 
matter fltoh| 1 970[ |Bodmer| 197 1 [ |Witten| 1 984} . The novel propo- 
sition was made that during the spin-down evolution of the neu- 
tron star a detonative (Niebergal et al. 2010) [Ouyed et aL]|201 1) 
phase transition to up-down-strange triplets would eject the outer 
layers of the neutron star at ultra-relativistic velocities (Keranen 
|et al.|2005||Ouyed & Leahy|2009| > (see the first panel of Fig. [2f. 
Follow-up studies of neutrino and photon emission processes dur- 
ing the QN JVogt et al.|2004||Ouyed et al.|2009b) l have shown that 
these outermost layers (10 -4 - 10 _3 M Q ) can be ejected with up to 
10 53 erg in kinetic energy. Nucleosynthesis simulations of the evo- 



lution of the neutron-rich QN ejecta (QNE) was found to produce 
primarily heavy elements with mass number, A > 1 30 ( Jaikumar et 
|al.|2007) . 

If the time (/qn) between SN and QN explosions is lengthy 
the SNE will have dissipated such that the QN essentially erupts 
in isolation. However, when ?qn is on the order of days a violent 
collision occurs reheating the extended SNE to a temperature of 
~ 10 9 K |Leahy & Ouyed|2008||Ouyed et al.|2009a| . The brilliant 
radiance of the re-shocked SNE (RSNE) fades as the photosphere 
recedes, eventually revealing a mixture of the inner RSNE and the 
deposited QNE material. In this work we revisit and extend that 
of |Leahy & Ouye d ( 2008 ) by introducing effects from the QNE to 
both the light curve and spectra of SN 2006gy. 



2.1 The Re-shocked Supernova Envelope 

The collision between the QNE and SNE creates a shock that prop- 
agates at a speed v shock through the SNE (see the middle panel of 
Fig. |2j. The shock breaks out of the SNE at a distance /?rsne,o at 
time 
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The RSNE which has a mass M RSNE is therefore fully re-shocked 
at a radius 

R * + v SN,max f QN 
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where R, is the radius of the progenitor star and v SN mlix is the max- 
imum velocity of the homologously (v oc r) expanding SNE. The 
outer edge of the RSNE grows as 



Rr 
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The core of the spherical RSNE will cool adiabatically from its 
initial shock temperature (T ) as, 
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For this proof-of-principle analysis we consider the temperature to 
follow a power law profile as: T RSNE (r, t) oc T^^ E (t) x r~P, where fj 
is a constant. 

To ensure that energy is conserved as continuum radiation is 
emitted by the hot RSNE (hRSNE) an equivalent amount of thermal 
energy is removed. Since the hRSNE is optically thick the radia- 
tive cooling process starts from the outer edge and a photosphere 
propagates inward leaving behind a cold outer layer (cRSNE), as 
shown in the right panel of Fig. [2] All species within the cRSNE 
are neutral and thus efficient absorbers of high energy radiation. 
While we leave a rigorous treatment of X-ray absorption for future 
work, modelling column density of the cRSNE (see Fig. |4]l yields 
an explanation for why CHANDRA observations show suppressed 
X-ray production by SN 2006gy l |Smith et al.|2007] >. 

2.2 The Hot Plate 

2.2.1 Formation 

The QN explosion sends ejecta propelling through space at speeds 
approaching c. The ejecta will eventually catch up to and collide 
with the SNE material and mix at the inner edge of the RSNE. The 
resulting composite shell will be referred to as the hot plate (HP) 
which has a constant thickness A^ H p- Conservation of momentum 
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tells us the HP will coast once it has captured all the SNE material 
it is capable of sweeping up. The final momentum of the HP will 
be equal to the momentum of the QNE plus that of all the SNE 
material (represented by a velocity space integral) swept into the 
HP: 
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where £qn is the energy of the QNE, Mhp is the mass of the 
HP and v H p is the final coasting speed of the HP. Assuming homolo- 
gous expansion (r oc v) and a constant density (p = M t /(4/37r^ lax )) 
the mass interior to velocity v is: 



Msne(v) = M t - 



(6) 



where M t is the total mass of the progenitor and vsN.max is the 
velocity at the outer edge of the SNE. Substituting equation [6] into 
[5]and solving for the velocity of the HP gives: 
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and using equation [7] in [6] assuming Msne(vhp) s> Mqne, we 
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Several factors such as gravitational braking, reverse shocks 
and mixing will lead to lower vhp values than that given by equation 
[7] A lower limit for vhp is given by the escape speed vnp.esc at the 
initial radius of the HP, #hp,o- The time, ? coast , for the HP to form 
and start coasting can be estimated as follows: 

The time taken for the QNE to catch up to and sweep up the 



HP material is (f coast - fQN 
be expressed as: 

^Hp(fcoast) = ttC(? coast - ?qn) 



), so the radius of the HP at time f coast can 
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where ac represents an average speed of the QN ejecta dur- 
ing the catch-up and sweeping phases. We explored a large range 
physical parameters of the QNE and found that typically the HP 
enters the coasting phase between 1.5 fQN and 3 ?qn, thus we chose 
as a fiducial value f coast ~ 2 fg N after which the radius of the HP is 
described by: 

^HP(0 = ^HP.O + VHP X (fQN - fcoast), 

2.2.2 Cooling 
For an adiabatic gas we have: 
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Since the HP is a coasting shell of constant thickness its vol- 
ume is given by 47rr 2 A,/?HP- Using y = | with equation[IT]thus gives 
the temperature as a function of time: 



T(t) = T 



Rw(t) 



(12) 



Photon diffusion will become important in the cooling of the 
HP when the HP and the overlying envelope become sufficiently 
transparent. We can approximate when this will happen (f crit ) by 
setting the mean free path to a fraction of the envelope size: 



A ~ 
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Up to a radius of q^sne we assume heating by the HP com- 
pletely ionizes the hydrogen in the envelope at all times f < f clit . 
Beyond this, the envelope may or may not be ionized depending on 
its temperature. The radius of the envelope at time f cr j, is given by 



^SNE( f crit) 



jfdit which leads to (using equation |13p 
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4000kms~ 1 ,a = 0.01,M e , 



25M 

and fi = 1.2 we get f crit as 200 days. Once f cril is reached, dif- 
fusion will quickly cool the HP within ~ 1 day. This calculation 
demonstrates that the HP is able to retain its heat for hundreds of 
days. Recognition of the sustained high temperature of the HP is 
essential to deciphering the enigma of SN 2006gy as this naturally 
describes both the late stage plateau in the light curve and the un- 
yielding broad component of the Ha line. 



2.2.3 Emission 

Radiation from the HP is treated as diffusion luminosity, given by: 
dD 



^^HP.p^Hpfv Ar : 



dt 



(15) 



where c v is the specific heat, ;i H p is the number density of the HP, 
AT represents the temperature difference at the interface between 
the HP and RSNE. The HP photosphere (Rhp, p ) defines the emitting 



surface, which moves at the diffusion rate, ^% (Leahy & Ouyed 



2008). 



2.3 Radiative Transfer 

The geometry of our model is quite simple as it contains only three 
components; a thin inner shell (the HP) surrounded by very thick 
outer shells (the hRSNE and cRSNE). For our analysis we use the 
astrophysical modelling software SHAPB^](Steffen et al.|2010 1 to 
build a geometric model of the RSNE and HP. Alongside being 
able to quickly and accurately construct 3-D geometries, density 
and temperature profiles, SHAPE is able to perform fast radiative 
transfer calculations. 

Radiative transfer in SHAPE is accomplished through a ray- 
casting algorithm. The aforementioned 3D model is subdivided 
into a 128 3 -cell grid, with each cell representing a volume (voxel) 
through which radiative transfer is performed. Given a temperature 
and density distribution, SHAPE calculates emission and absorp- 
tion coefficients for each voxel in the grid. A ray is then cast from 
behind the grid to the observer, performing radiative transfer cal- 
culations along each voxel for each frequency band. The emerging 
ray can then be used in the generation of the light curve and spec- 
trum plots. This allows us to simultaneously produce light curves 
and spectra which can then be compared to the observations of SN 
2006gy. For this analysis we use rays consisting of 150 bands rang- 
ing from 3 X 10 14 Hz - 3 x 10 18 Hz. The plot of the light curve only 
uses R-band frequencies from ~ 3.76 x 10 14 Hz - 5.76 x 10 14 Hz, 
but the energy conservation requires the entire range. For the spec- 
tra, we use 200 bands from 4.4118 x 10 14 Hz - 4.6875 x 10 l4 Hz 
(spanning the frequency of Ho-). 
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Unlike most radiative transfer codes, SHAPE does not calcu- 
late the temperature of each cell in the grid. Instead, the tempera- 
ture is an input in the form of an analytic distribution. Although less 
accurate than a full radiative transfer treatment, the speed gain al- 
lows for an interactive parameter search as required by this work. In 
summary, SHAPE answers the question "Given a temperature and 
density, what is the intensity of the light emitted by an object?". 



2.4 Model Parameters 

The parameters used to construct the light curve and spectra of 
SN 2006gy can be divided into three categories: Constrained, Free 
Physical, and Free Radiative Transfer. 



2.4. 1 Constrained Parameters 

The following parameters that are used in our model are con- 
strained by observations. 

• M t : the total mass (M RS ne + ^hp) is constrained from above 
by ~ 40M G (anything heavier will likely form a black hole rather 
than a neutron star ( |Heger et al. {2003) ) and ~ 20M o from below 
(anything less massive will not form a quark star ( |Staff et aL|2 006 ))■ 

• VsN,max : the outer edge velocity of the expanding RSNE is con- 
strained to be 4000 km s 1 by the blue-side absorption feature (as 
noted in Fig.|9| observed in the Ho- spectrum (Smith et al.|2010| >. 

• T : the initial temperature of the RSNE and HP is constrained 
by fitting the broad component of the Ho- line for each spectrum 
assuming it is thermally broadened. The extrapolated To value is 
on the order of ~ lxl0 9 K. 

• v s hock : the speed of the shock as it travels through the RSNE 
is ~6000 km-s~', this is found using shock physics as kTo = 



at an emission mechanism. In addition to the absorption in equation 
|17| we include a Thomson scattering term of the form: 



" ™ 1,2 where the initial temperature (T ) is known 



2.4.2 Free Physical Parameters 

In our model the SN explosion of a massive star is followed a time 
?on later by a QN explosion. The QNE/SNE collision sends a shock 
through the SNE creating two distinct layers: the HP and RSNE. 
The HP has a mass Mhp, coasting speed vhp, initial radius /?hp.o 
and constant thickness A_R H p- The RSNE has a mass, M RS ne an d 
temperature profile defined by the power law index, B. We assume 
constant density for all components. These variables make up the 
free physical parameters in our model which we find by fitting our 
theoretical light curve to the observations of SN 2006gy. 



2.4.3 Free Radiative Transfer Parameters 

The radiative transfer in SHAPE requires the specification of an ab- 
sorption and emission coefficient for each frequency, v. The emis- 
sion coefficient used in this analysis is of the form: 



An] 



.liv/kT 



(16) 



Where A is a multiplicative factor. The corresponding absorp- 
tion coefficient (assuming LTE) is: 



h_ 

By 



(17) 



K V ,TH = Bn.Cr^H 



(18) 



Where <r TH is the Thomson scattering cross section and B is a 
multiplicative factor. Due to the computational and time demands 
of treating multiple scattering, we simply assume that scattering is 
isotropic and everything removed from the beam will eventually 
be put back in. Of course the scattering is not perfectly elastic and 
energy will be lost. Therefore B represents the fraction of scattered 
light that will NOT be scattered back into the beam. 



3 RESULTS 
3.1 Photometry 

Through the exploration of parameter space using SHAPE, we find 
that the light curve is well reproduced by a progenitor mass ranging 
from 20-40 M a . We have fit three models with progenitor masses 
of 20, 30 and 40 M Q with the parameters of each fit listed in table 
[T] Each uses an envelope temperature profile of B = 0.2 with an 
initial temperature of T = 0.9 x 10 9 K. The initial energy of the QN 
ejecta is set to ~ 1 x 10 52 ergs ( |Keranen et al.|2005 I We find that 
the actual velocity of the HP differs from the theoretical (equation 
|7]( by a factor of 0.1. As previously mentioned, factors such as 
gravitational braking, reverse shocks and mixing will lead to lower 
Vhp values. 

The emission is simulated using the coefficients of equations 
[16] [17] and [18] with the A values listed in table [T] and B = 5 x 
10~ 4 . While the shock is still moving through the SNE (t < f SB0 ) 
detailed shock physics would be required to properly describe the 
corresponding light curve, thus we begin our photometric analysis 
at time t = t SBO ^ 30 days. Our SN explosion date is 5, 10 and 20 
days earlier than that inferred by Smith et al. (2010) for the 40, 30 
and 20 M B fits respectively. 

Figure [3] shows the fits for the various progenitor masses, all 
of which match the data well. Figure [4] shows how the different 
components contribute to the light curve. The high initial tempera- 
ture and radius of the RSNE provide ample luminosity (see the blue 
dash-dot line in Fig. EJ to explain the broad peak of the observed 
light curve. As the photosphere recedes, radiation from the under- 
lying HP becomes more prominent and the drop in the overall light 
curve (red solid line in Fig. HI is mitigated by support from the HP 
(green dashed line in Fig.Rj). The slowly coasting and thus slowly 
cooling HP provides a natural explanation for the plateau observed 
in the late stage of the SN 2006gy light curve. 



Where B v is the Planck function. The form of equation [T6] is 
similar to that of bound-free continuum emission, perhaps hinting 



3.2 Spectroscopy 

For our proof-of-principle spectroscopic analysis we aimed to de- 
termine the cause of the overall structure of the Ha spectrum of 
SN 2006gy, as to date only a phenomenological analysis has been 
done (ie. fitting a Gaussian versus a Lorentzian profile ( Smi th et aT] 
|2010| [Chatzopou los et al.|2011[ l). Displayed in Fig.[6]is our model 
Ha line (red solid line) compared to spectral observations of SN 
2006gy for days 36, 65, 71, 96, 125, 154, 179 and 209. As shown 
in the right panel of Fig. [2] after shock break out our model con- 
tains three components; an outer layer of cold RSNE (cRSNE), a 
hot region of RSNE (hRSNE) and an innermost thin HP shell. Each 
of these constituents contribute to the overall shape of the Ha line. 
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For our modelling we ignore contributions from the two identified 
[N II] lines at 6548 and 6583 Angstroms. 

By viewing the width of the broad Ho- component in the light 
of thermal rather than Doppler broadening, the reason for the curi- 
ous late stage breadth becomes clear. Considering the slowly coast- 
ing HP as the thermal source, the width of the broad component is 
then proportional to Tap. Plotted with a green dash-dot line in Fig. 
[9]is the broad Ho- component for day 96 (width oc r HP = 7.5 x 10 8 
K). We determined the temperature evolution of the HP by fitting a 
Gaussian to the broad component of the Ho- line. The decrease in 
temperature agrees remarkably well with the predicted HP temper- 
ature from our light curve (top-right panel of Fig. [5}. 

The peculiar asymmetry of the Ho- line observed in SN 2006gy 
is inescapable in our scenario. The wide blue-side absorption fea- 
ture that does not dip below the continuum is due to RSNE material 
along our line of sight absorbing emission from the underlying HP. 
The blue-side feature extends out 4000 km-s~' as both the hRSNE 
and cRSNE contribute to the absorption. The blue-side absorption 
is strongest at times prior to day 96 and then becomes weaker as 
the absorbing RSNE material becomes increasingly diffuse and in- 
creasingly cool. The bottom right panel of Fig.[6]most clearly dis- 
plays the evolution of the blue-side absorption in our model. 

At early times the narrow Ho- peak is broadened by thermal 
motion. At some time between day 71 and 96 the temperature of 
the hRSNE drops such that kinematic Doppler broadening takes 
over. As the location of the Ho- emission region recedes into slower 
moving RSNE material, emission from the hRSNE originates from 
ever lower velocities. This accounts for the discontinuity in the Ho- 
line at -1200 km-s~' on day 96 seen in Fig. [9] However the cold 
outer layer of the RSNE continues to absorb out to 4000 km-s~' 
(see green dashed line in Fig. [9}. 

The sharp P-Cygni profile located at the center of the Ha line 
seen in observations on days 96, 125 and 154 is characteristic of an 
expanding (velocity « 120 km-s~') shell of Hydrogen. This veloc- 
ity suggests the sharp P-Cygni feature is likely due to stellar wind 
from the progenitor star. This feature does not seem to be typical 
for super-luminous SNe as it has not been observed in the super- 
luminous SNe: SN 2008am, ( |Chatzopoulos et al.|201l| >, SN 2008fz 
( [Drake et al.|[20T0> , SN 2003ma ( |Rest et al.poil^ SN 2008es 
( |Miller et al.||2009l |Gezari et al.|[2009l > or SN 2007bi jGal-Yam| 
|et al.|2009| >. 

For illustrative purposes we added an outermost shell of 
Hydrogen to our model in order to simulate a stellar wind 
(p v r 2 =const.) from the progenitor. Our low mass (« 0.5 M G ) stel- 
lar wind extends from the outer edge of the cRSNE to a radius of 
3 x 10 14 m and has a temperature of 8000 K. Fig. [7] displays a 
comparison of the observed evolution of the Ho- line of SN 2006gy 
with our (stellar wind included) QN model. With the exception of 
the stellar wind, the same model parameters were used to create the 
fits seen in Figures|6]and|7] As seen in Fig.|7]at early stages (days 
36, 65 and 71) the stellar wind continuum is too low for absorption 
to occur and the stellar wind adds a narrow peak to the Ho- line. The 
continuum emission of the stellar wind then increases to a level that 
a strong P-Cygni profile becomes present by day 96. The P-Cygni 
profile diminishes with time as the strength of the stellar wind con- 
tinuum falls. By day 209 the P-Cygni profile is gone leaving only 
a narrow emission peak as the stellar wind contribution to the Ho- 
line. The bulk structure of the Ha line can be described by emis- 
sion from the HP and the RSNE (see fig.[6](, while the surrounding 
low mass stellar wind adds the fine details to the spectral line (see 
Fig.|7J. A side by side comparison of our model Ho line with the 
observed spectral line can be seen in Fig. [8] 



In this proof-of-principle analysis we found that the Ho- emis- 
sion originated from a region interior to the photosphere (see Sect. 
2.1 for definition of photosphere) at all times in the evolution of the 
spectrum. This is to be expected as the photosphere defines the lo- 
cation of continuum emission while Ho- emission occurs at a chro- 
mospheric region which is dependent on the density and tempera- 
ture of the plasma (Mihalas 19781. As future work, improvements 
must be made to our formalism to include such elements as asym- 
metries since higher order effects could have an impact on the spec- 
troscopic behaviour of this system. A more detailed study of the 
structure of the RSNE and the role of the density and temperature 
profiles on the spectroscopy of the QN will as well be undertaken 
in future studies. 



4 DISCUSSION 

The HP temperature (Tup) provides us with an opportunity to 
check our model for consistency. Comparing our model 7hpW = 
T ' (red solid line in the top-right panel of Fig. jij with 

values found by fits to the Ho- line (black circles plotted in top-right 
panel of Fig. [5} we can see that they are in good agreement. 

As noted by Smith et al. (20101 there exists a second peak of 
Ho- luminosity accompanied by peaks in Fe II and Ca II luminos- 
ity that evolve in concert (see the bottom panel of Fig. |5j. This is 
expected in our QN scenario, since during the creation of the HP, 
ploughing of the inner RSNE material will mix the elements lead- 
ing to their simultaneous evolution. As seen in the bottom panel 
of Fig. [5] the rise of the Fe II (red diamonds), Ca II (blue squares) 
and Ha luminosity (black circles) correspond with the emergence 
of the HP (green dashed line). 

During the formation of the HP, elements with A > 130 (in the 
QNE, Jaikumar et al. ( 2007)) will mix with the lighter elements in 
the SNE. We predict some A > 130 elements should appear late in 
the evolution of the SN 2006gy spectrum evolving in concert with 
the Fe II and Ca II. The intensity of the A > 130 emission lines 
should be much lower than those from the RSNE as 



M, 



QNE 



M H 



< 10- 



(19) 



Plotted as a solid black line in Fig.[4]is the column density of 
the cRSNE. After shock breakout the column density of the neutral 
cRSNE rapidly rises to 10 25 cnT 2 , thus quickly becoming opaque to 
X-rays. This naturally explains why CHANDRA observations show 
minimal X-ray production from SN 2006gy. 

We were able to reproduce the light curve of SN 2006gy us- 
ing a range of progenitor masses. The lower mass models required 
a higher emission coefficient (A value) and longer time delays to 
recover the fit. Because T does not vary from one fit to another 
and the continuum levels remain the same (hence the same fit to 
the light curve), the spectra are unaffected by varying the above pa- 
rameters. The range of possible A coefficients highlights our igno- 
rance of the emission process and more detailed radiative transfer 
calculations will be required to narrow down the progenitor mass 
and the time delay between SN and QN. Regardless, our models 
do show that given a progenitor mass between 20M o and 40M o , 
enough energy is available to explain the amplitude and breadth of 
the SN 2006gy light curve. 
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5 CONCLUSION 

The QN model presented in this proof-of-principle analysis simul- 
taneously provides a very good fit to both the light curve and spec- 
tra of SN 2006gy. The broad, luminous peak of the light curve is 
described by radiation emitted from the hRSNE. The late stage 
plateau in the light curve, which is thus far inexplicable in other 
models, is naturally explained by emission from the slowly coast- 
ing HP. These distinctive features of the QN model also account 
for the peculiar structure and evolution of the Ha line. The persis- 
tent broad Ho- component, which is not understood in the context of 
any other model, is uniquely described by emission from the slowly 
cooling HP of the QN model. The blue-side absorption feature ob- 
served in the Ha line is caused by both the cRSNE and hRSNE 
which attenuate emission from the underlying HP. As the hRSNE 
radiates energy it builds a cold outer layer (the cRSNE) which is 
opaque to high energy radiation. This illustrates how an explosion 
as energetic as SN 2006gy can be surprisingly quiet in X-rays. The 
QN is currently being used to model other super-luminous SNe. 
The findings will be published in an upcoming paper. 



6 FUTURE WORK 

This work is intended to demonstrate that a QN explosion of a neu- 
tron star only days after the initial SN is able to account for the 
light curve and spectra of SN 2006gy. However, the exact emission 
mechanism of the radiation is not well understood. The emission 
coefficient that best reproduces the rise and fall of the light curve 
resembles bound-free continuum emission (boosted by a factor of 
A). The boost may be due to the contribution of some other emis- 
sion mechanism such as free-free. Future work will concentrate on 
a more detailed treatment of radiative transfer than possible with 
the current computer code. 
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Table 1. Model parameters for different progenitor masses (M t ). Parameters listed with * are determined through equations and are not adjustable. 



M,(M G ) 


M env (M )* 


M nv (M e y 


fQN (days) 


A 


Rhp.o (cm)* 


ARhp (cm) 


v H p (kms ')* 


40 


36.3 


3.7 


10.5 


1750 


3.3 x 10 14 


2.5 X 10 12 


1809.7 


30 


26.6 


3.4 


11.5 


4000 


3.9 x 10 14 


5.0 X 10 12 


1944.7 


20 


16.9 


3.1 


12.5 


12000 


4.6 x 10 14 


10.0 X 1() 12 


2152.1 




Figure 1. Photometric and spectroscopic observations of SN 2006gy. Left: The light curve (R-band) of SN 2006gy (black squares) is compared to light 
curves of three supernovae; a typical type la (SN 1998dh, blue long dash line), a typical type II (SN 1999em, green short dash line) and SN 1987a (mauve 
dash-dot line). SN 2006gy is ~ 10 times brighter than a typical type la and ~ 100 times brighter than a typical type II. As well, the light curve displays an 
unusual plateau beginning at ~ day 200. Right: Spectral observations of the SN 2006gy Ha line for days 96 (black solid line), 125 (blue long dash line), 154 
(red short dash line) and 179 (green dash-dot line) after inferred explosion date are overplotted and normalized to the same continuum level. In this figure and 
all subsequent spectroscopic figures, the inferred explosion date is that cited by Smith et al. 2010. The Ha line displays a persistent broad component and a 
blue-side absorption feature that decreases with time. 
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Figure 2. Schematic representation of the quark nova scenario. Displayed are cross sections of our model at three different phases. Left: The instant the 
QN detonates inside the expanding envelope of SNE. At this time the QNE is expelled at ultra-relativistic velocities leaving behind a quark star (QS). Middle: 
The interval in which the SNE is being shocked by the collision with the QNE. Meanwhile at the inner edge of the envelope RSNE and QNE material are 
mixing to form the HP. Right: The evolution of the QN scenario after shock breakout. Radiative cooling builds a cold outer layer of the RSNE (cRSNE) as the 
photosphere recedes inward. At this time the HP is fully formed and coasting inside the RSNE. 
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Figure 3. Comparison of models with different progenitor mass. The blue dash-dot line shows the R-band light curve fit for M, = 40M o , the red solid 
line for M, = 30M Q , and the green dashed line for M, = 20M o . The parameters for each curve are given in table[T] The open black squares represent R-band 
observations of SN 2006gy over 230 days beginning at the inferred explosion date of Smith et al. (20071. The different starting times for each curve reflect the 
different shock breakout times which is later for longer time delays. 
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Figure 4. Comparison of the QN model to the R-band light curve of SN 2006gy. The red solid line shows the 30 M Q best fit R-band light curve for our QN 
scenario which includes the RSNE and the HP. Also plotted are the individual RSNE (blue dash-dot line) and HP (green dashed line) component R-band light 
curves as well as the column density of the cold outer layer of the RSNE (cRSNE) (black solid line). The open black squares represent the R-band observations 
of SN 2006gy over 230 days beginning at the inferred explosion date of Smith et al. ( 2007). 
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Figure 5. Further signatures of the HP. Left: In this figure our predicted temperature of the HP is compared with values determined from fits to spectral 
observations of SN 2006gy. The red solid line represents the temperature of the HP from our best fit light curve (30M o ). Plotted as black closed circles are 
the HP temperatures found by fitting a Gaussian to the broad component of the Ho- line. The error bars represent a 5 % measurement error. Right: This figure 
shows the evolution of the luminosities for the Ha line, Ca II line at 8662 Angstroms and the red Fe II blend over roughly 7 1 30-7628 Angstroms from Smith 
et al 2010. After day 90 the Fe II and Ca II emission track the sudden peak in Ha luminosity almost exactly. This corresponds to the time when the opacity of 
the RSNE has dropped enough such that HP emission becomes a strong contributor to the overall light curve of SN 2006gy. For comparison the R-band HP 
luminosity from our model light curve is plotted as a green dotted line. 
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Figure 6. Evolution of the Ha line in SN 2006gy. This series displays the evolution of the observed Ha line (grey dashed line) from SN 20()6gy (data from 
Smith et al. 2010) compared with our model Ha line (red solid line). The dates shown are: 36, 65, 71, 96, 125, 154, 179 and 209 days after the inferred 
explosion date (left to right, top to bottom). The narrow peak in our model Ha line is due to emission from the hRSNE. The blue-side absorption feature is 
caused by the entire RSNE which has an outer edge velocity of 4000 km s . The broad component of the Ha line is emission from the HP which is coasting 
at a velocity of 95 km s and thus slowly cooling from its initial shock temperature (0.9 X 10 9 K). As the Ha line evolves, the component due to the HP 
becomes more prominent and the blue-side absorption feature diminishes. Displayed in the bottom right panel is our model Ha line produced on days: 96, 
125, 154 and 179 normalized to the same continuum level and overplotted. The blue-side of the broad component shows strong absorption that decreases with 
time due to decreasing RSNE temperature and density, while the red-side remains notably constant. 
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Figure 7. Evolution of the Ha line in SN 2006gy including stellar wind. This series displays the evolution of the observed Ha line (grey dashed line) from 
SN 2006gy (data from Smith et al. 2010) compared with the Ha line (red solid line) from our model which includes stellar wind from the progenitor. The dates 
shown are: 36, 65, 71, 96, 125, 154, 179 and 209 days after the inferred explosion date (left to right, top to bottom). The intensity of the continuum emission 
from the wind is initially too low to allow for absorption. During intermediate times the stellar wind continuum has increased such that a P-Cygni profile 
occurs, which is most prominent at day 96 and diminishes there after. At late stages the stellar wind continuum is once again too low to allow for absorption. 
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Figure 8. Ha line over-plotted. For both panels of this figure the Ha line from days 96, 125, 154 and 179 are normalized to the same continuum level and 
over-plotted. Left: Displays the spectral line produced by our model. Right: Displays the observations of SN 2006gy (data from Smith et al. 2010). 
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Figure 9. Features of the Ha line. This figure shows the observed Ho- spectrum for day 96 (grey solid line) compared with the unabsorbed broad component 
from the HP (blue dash-dotted line) as well as our model Ho- line (red solid line) and a version that does not include the cold outer layer of the RSNE (cRSNE) 
(green dashed line). The velocity of the outer edge of the RSNE and the RSNE photosphere are indicated with arrows. This figure clearly demonstrates the 
contribution of our three components (HP, hRSNE and cRSNE) to the structure of the Ha line. 



